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1. Introduction 

The principle goal of the relativistic heavy-ion collision experiments is the discovery of 
new state of matter, the so called quark gluon plasma (QGP). The various measurements 
taken at CERN SPS and at BNL RHIC do lead to strong 'circumstantial evidence' for the 
formation of the QGP [1,2]. Evidence is circumstantial as any direct formation of the QGP 
cannot be identified. Only by some noble indirect diagnostic probes like the suppression of 
the J /if! particle, the jet quenching, the enhanced production of strange particles, specially 
strange antibaryons, excess production of photons and dileptons, etc. the discovery can be 
achieved. 

In order to understand the properties of a QGP and to make unambiguous predic- 
tions about signature of QGP formation, one needs a profound description of QGP. For 
this purpose we have to use QCD at finite temperature and chemical potential. There 
are two different approaches: 1) Lattice QCD is a non-perturbative method for solv- 
ing the QCD equations numerically on a 4-dimensional space-time lattice [3]. In this 
way all temperatures from below to above the phase transition are accessible. 2) Pertur- 
bative QCD at finite temperature [4,5] is based on the fact that the temperature depen- 
dent running coupling constant is small at high temperatures due to asymptotic freedom, 
T — > oo =>■ a s (T) = g 2 /Air — > 0. At a typical temperature of T = 250 MeV we expect 
a s = 0.3 - 0.5. This suggests that perturbation theory could work at least qualitatively. 
This corresponds to an expansion in a s , which can be performed conveniently by using 
Feynman diagrams for partonic scatterings. 

2. Thermal Field Theory 

Thermal field theory is a combination of three basic branches of modern physics, viz., 
quantum mechanics, theory of relativity and statistical physics. Our aim is to derive Feyn- 



1 



M. G. Mustafa 



man diagrams and rules at T > (p 7^ 0). We start by repeating some of the basic facts 
in equilibrium statistical mechanics. The statistical behaviour of a quantum system, in 
thermal equilibrium, is studied through an appropriate ensemble and is defined as 

p{P) = e~ m , (1) 

where f3 = 1/T (the Boltzmann constant ks = 1) and H is appropriate Hamiltonian for 
a given choice of ensemble. Given the density matrix, the finite temperature behaviour of 
any theory is specified by the partition function 

Z((3) = Trp(/3) = Tre-f™ . (2) 

The thermal expectation value of any physical observable can be written as 

(A) (j = Z- 1 (f3)TT[p((3)A]=Z-\f3)TT[e-^A\ , (3) 

and that of correlation function of any two observables is given as 

(AB)/} = Z~\P)Tr [ P {p)AB] = Z" 1 (/?)■& [e^ n AB] , (4) 

For a given Schrodinger operator, A, the Heisenberg operator, Ah {t) can be written as 

A H (t) = e mt Ae- mt . (5) 

The thermal correlation function of two operators can also be written as 

{A H {t)B H {t'))p = Z-\I3)Tt [e- 0H A H (t)B H (t')] - (B H (t')A H (t + 0))p , (6) 

which holds irrespective of Grassmann parities of the operators. Eq.(6) is known as the 
Kubo-Martin-Schwinger (KMS) relations and will lead to (anti)periodicity in various 2- 
points functions at finite temperature. 

We would like to note that in (2) "Tr" indicates the sum over expectation values in all 
possible states in the Hilbert space and there is an infinite number of such basis in quantum 
field theory [6,7]. The partition function of a statistical system cannot be computed ex- 
actly even if one makes a perturbative expansion into a power series in coupling constant, 
g. The Matsubara formalism, known as also imaginary time formalism, provides a dia- 
grammatic way of computing the partition function and other relevant physical quantities 
perturbatively. 

2.1 Matsubara formalism 

The Hamiltonian of a system can be decomposed as 

H = H + H' , (7) 

where Ho and H' are the free and interaction parts, respectively. The density matrix in (1) 
becomes 

p{p) = po(/W) with p (p) = ■ S(J3) = e^ n . (8) 
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The density matrix can have evolution equation with < r < j3: 

= -Hopo(r) and ^ = -Hp(r) = - (Ho + W) p(r) . (9) 
Now S(t) satisfies the evolution equation, following (8) and (9), as 

= dT P( T ) + Po -q^- = -W / (r)5(r) , (10) 

with H'j(t) = pQ 1 (t)H' po(r) — e Tno H'e~ Tn " in a modified interaction picture which 
is similar to zero temperature field theory. Note that for a real r such transformation may 
not be necessarily unitary but for imaginary values of r = it it will be so. This makes 
Matsubara formalism a imaginary time formalism. It is also to be noted that under such 
rotation to imaginary time the field remains hermitian with the appropriate definition of 
hermiticity for complex coordinates [6], viz., ${z) — 4>{z*). Integrating (10) one can 
obtain [6] 

S(fi)=T T (V So ' dTU '>^ , (11) 

where T T is the ordering of r variable. The following points to be noted: 1) Eq.(l 1) resem- 
bles the 5-matrix in zero temperature field theory with the exception that the time integra- 
tion becomes now finite. 2) Now one can expand the exponential in (1 1) into a power series 
of the coupling constant and each term would then correspond to a Feynman diagram. 3) 
In this formalism Wick's theorem can easily be generalized to finite temperature. 



2.2 Matsubara frequency 

The 2-point Greens functions can be defined in Heisenberg picture as 

^(ry) = (T T (^(r)4(r')))^ = 2;- 1 ( / 9)'&e-' 3w r T (^(r)4(r')) , (12) 
with (f> is a complex field and the r ordering is defined as 

T T (<Mt)4(t')) = 0(r - r')^(r)4(r') T 0(r' - t)^ h {t')4> h {t) , (13) 

where negative sign in the second term is for fermionic field as the t ordering is sensitive 
to the Grassmann parity of the fields. The spatial dependence and spinorial indices are 
omitted for the time being as they are similar to the zero temperature field theory. 

Exploiting the cyclicity properties of the trace and (13), one can show that for r > 
(12) satisfies (anti)periodic conditions 

G(0, r) = T = TGpifi, t) , (14) 

for (fermion)boson which restrict the time in finite interval [0, 0\. The (14) is equivalent to 
KMS relation defined in (6) with a suitable rotation to imaginary times. 
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The Fourier transform of 2-point functions would involve the discrete frequencies as 
those are defined in finite time interval and in general can be written as 



^^Ee-^^K) , d5) 

n 

1 r 

QpM = 2 J dre^ r g (T) , (16) 



where u) n = mr/(3 with n = 0, ±1, ±2, • • •. The (anti)periodicity conditions in (14) 
restricts (16) as 



-0 

= \ (1 T (-1)") J Are^g p (r) , (17) 

...... r . f 2mrT for bosons . . , 

yielding discrete frequencies u) n — < _|_ for fermions ' known as Matsub- 

ara frequencies. 

Two important consequences of imaginary time formalism: 1) The time r is restricted 
to the finite interval [0, (3] due to (anti)periodic conditions for (fermion)boson. 2) The 
finite time interval changes the integral over the time component of momentum at zero 
temperature to a Fourier series over Matsubara frequencies. 



2.3 Propagator for any theory 

Since the spacial dependencies of 2-point functions are continuous like zero temperature 
field theory, therefore, putting all the coordinates together one can write 

SfMr) = ^E/ ?|^ e ~ i(W " T ~^(^n) , (18) 

g p (k,u n ) = J dr J dy^-^^fiV) . (19) 

The zero temperature Greens function satisfies the bosonic Klein Gordon equation as 

(d^ + m 2 ))g(x) =5 i (x), (20) 

with the choice of metric in Minkowski space is (+, — , — , — ). While rotating to the imag- 
inary time (t — > — it => g — ► —g$) the finite temperature Greens function satisfies the 
equation 

(jj + y2 - ™ 2 ) T ) = S 3 (x)6(r) ■ (21) 
On substituting (18) in (21), the momentum space propagator is obtained as 
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uj 2 + k 2 + m 2 



(22) 



Now substituting (22) in (18) and performing the frequency sum [5], one can obtain 



where uik = v'k 2 + m 2 and nsi^k) = (e ,3 " fc — is distribution function for boson. 
Now putting t = it and nsi^k) = 0, one can recover the zero temperature 2-point Greens 
function for this case. Similarly, one can also obtain 2-point Greens function for fermion 
from Dirac equation. 

2.4 Feynman rules for finite temperature field theory 

1 . Replace T = propagator by T ^ propagator as obtained which carry T depen- 
dence through Matsubara frequency. 

2. The vertex is same as T = case. 

3. Replace loop integrals by / d A K/ (27r) 4 — > T J2 n J d 3 k/(2ir) 3 . 

4. Symmetry factors are same as T = case. 

3. Hard Thermal Loop (HTL) Approximation 

Restricting only to bare propagators and vertices perturbative QCD can lead to serious 
problems, i.e., infrared divergent and gauge dependent results for physical quantities, e.g., 
in the case of the damping rate of a long wave, collective gluon mode in the QGR The sign 
and magnitude of the gluon damping rate was found to be strongly gauge dependent [8,9]. 
The reason for such behaviour is the fact that bare perturbative QCD at finite temperature is 
incomplete, i.e., higher order diagrams missing in bare perturbation theory can contribute 
to lower order in coupling constant. In order to overcome these problems, the HTL resum- 
mation technique, an improved perturbation theory, has been suggested by Braaten and 
Pisarski [10] and also by Frenkel and Taylor [11], in which those diagrams can be taken 
into account by resummation. The starting point is the separation of scales in the weak 
coupling limit, g « 1, since there are two momentum scales in a plasma of massless par- 
ticles: (i) hard, where the momentum <~ temperature, T and (ii) soft, where the momentum 
<~ thermal mass (~ gT, g « 1). 

3.1 HTL perturbation theory (HTLpt) 

The HTL resummation technique allows a systematic gauge invariant treatment of gauge 
theories at finite temperature and chemical potential taking into account medium effects 
such as Debye screening, effective quark masses, and Landau damping [10]. The gener- 
ating functional which generates the HTL Green functions between a quark pair and any 
number of gauge bosons can be written [10,12] as 




(23) 
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5£ = m ^ (k^d) ^ ' (24) 

where is a light like four-vector, D M is the covariant derivative, m q is thermal quark 
mass, and ( ) is the average over all possible directions over loop momenta. This functional 
is gauge symmetric and nonlocal and leads to the following Dirac equation 

pip = SV» + T^t/j + AT^A" V + • , (25) 

where we have suppressed the color index. In the HTL approximation the 2-point function, 
£ ~ gT, (quark self-energy) is of the same order as the tree level one, Sq 1 ^) ~ If ^ gT 
(in the weak coupling limit g <C 1), if the external momenta are soft, i.e. of the or- 
der gT. The 3-point function, i.e., the effective quark-gauge boson vertex, is given by 
gTfj, = g("in + STfj,), where ST^ is the HTL correction. The 4-point function, g 2 T^ v , 
does not exist at the tree level and only appears within the HTL approximation [10]. These 
Appoint functions, which are complicated functions of momenta and energies, are inter- 
related by Ward identities. The HTLpt [13] is an extension of the screened perturbation 
theory [14] in such a way that it amounts to a reorganization of the usual perturbation 
theory by adding and subtracting the nonlocal HTL term (24) in the action [10,11]. The 
added term together with the original QCD action is treated nonperturbatively as zeroth 
order whereas the subtracted one as perturbation. Equivalently, for calculating physical 
quantities this amounts to replace the bare propagators and vertices by resummed HTL 
Green functions [10]. However, HTLpt suffers from the fact that one uses HTL resummed 
Green functions also for hard momenta of the order T. A systematic description of phys- 
ical quantities requires an explicit separation of hard (~ T) and soft (~ gT) scales [10]. 
Following this approach, a number of relevant physical quantities, such as signatures of 
the QGP, has consistently been calculated [7,15-19]. 



3.2 Quark self energy in HTL- approximation 

— > — L — > — ^ — 

Figure 1. Quark self energy having quark momentum K. 

The most general ansatz for the fermionic self energy in rest frame of the plasma is given 
by [20] 

Z(K) = -a(k ,k)#-b(k ,k) 7o , (26) 

where K — (ko, k), k = |k| and the quark mass is neglected assuming that the temperature 
is much larger than the quark mass, which holds at least for u and d quarks. The scalar 
quantities a and b are given by the traces over self energies within HTL approximations in 
figure 1 as, respectively, 

a(k , k) = JLj [tr m) - k tr ( 7o S)] = ^ (l - | In ^±|) , (27) 
>(*>.*) = -h tr - k tr (#E)] = J [-ko + In g±| j , (28) 
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where m 2 q = g 2 T 2 /6, is the effective quark mass. The quark self energy in (26) has an 
imaginary part below the light cone, kg — k 2 < 0, representing Landau damping (LD) for 
space like quark momenta. Furthermore, the general ansatz in (26) is also chirally invariant 
in spite of the appearance of an effective quark mass. 



3.3 Effective propagators and vertices in HTL-approximation 



-> — + — > — L — => 



Figure 2. Effective quark propagator in HTL-approximation 

Resumming the quark self-energy by using the Dyson-Schwinger equation, the effective 
quark propagator in figure 2 can be written as 



-i _ 7o-k'7 . 7o + k-7 



2D+(K) 2D_(K) ' 
for massless quarks with a decomposition into helicity eigenstates, where 

D ± (k , k) = (1 + a(ko, *)) (-*o ± k) + b(k , k) . 



(29) 



(30) 




Figure 3. Quark dispersion relation in HTL-approximation along with the free 
(uj — k) one. 



The relevant QED like HTL-vertex is related to this propagator through the Ward iden- 
tity [19,11] 



K^{K 1 ,K 2 -K) = S- 1 {K 1 )-S~\K 2 ) 



(31) 



Now, the zeros of D±(K) give the in-medium propagation or quasiparticle (QP) dis- 
persion relation. As shown in figure 3 the upper curve 0J+(k) corresponds to the solution 
of D + {K) = 0, whereas the lower curve tJ-(k) represents the solution of D-(K) = 0. 
Both branches start from a common effective mass, m q , obtained in the k — > limit [19]. 
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The lu + (k) branch describes the propagation of an ordinary quark with thermal mass, and 
the ratio of its chirality to helicity is +1. On the other hand, the oj_ (k) branch corresponds 
to the propagation of a quark mode with a negative chirality to helicity ratio. This branch 
represents the plasmino mode which is absent in the vacuum but appears as consequence of 
the medium due to the broken Lorentz invariance, and has a shallow minimum. This cor- 
responds to a purely collective long wave-length mode, whose spectral strength decreases 
exponentially at high momenta. For high momenta, however, both branches approach the 
free dispersion relation. 

The spectral representation [7,19] of HTL propagator in (29) can now be written as 

p± = -Imf±(x + ie) = Polc(fc^ > k 2 ) + Cut(fc^ < k 2 ) . (32) 



4. Thermal Hadronic Correlation Function 



While lattice calculations of hadron properties in vacuum have reached quite satisfactory 
precision, little is known from such first principle calculations about basic hadronic prop- 
erties in a thermal medium. We want to analyse [21] its behaviour in the high temperature 
limit through the thermal correlation functions in terms of quark fields. 



4.1 Definition 

Meson correlators are constructed from meson currents Jm(t,x) = q(T,x)TMq(T 7 x), 
where Tm — 1. 75, 7/*, 7^75 for scalar, pseudo-scalar, vector and pseudo-vector channels, 
respectively. The thermal two-point functions in coordinate space, Gm(t, x), are defined 

as 

G m (t,x) = (J M (r,x)jl(0,0)) = T £ J e -i(«»T-jS0 XM ( Wn ,0) , (33) 

n— — oo 

where r <G [0, 1/T], and the Fourier transformed correlation function Xui^mP) is given 
at the discrete Matsubara modes, u„ = 2nirT. The imaginary part of the momentum space 
correlator gives the spectral function cjm{w,p), 

XM{uJ n ,p) = - ( duj °" M (^'^ . aM (<jj p) = llm xm(u ,p). (34) 

Using (33) and (34) we obtain the spectral representation of the thermal correlation func- 
tion in coordinate space at fixed momentum (J3 = 1/T), 

^ , f°° , , - cosh(w(r-/3/2)) 

Gm(t,p)= duj<T M (u;,fi \) 7 ;; ■ (35) 

Jo sinh(w/3/2) 
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Q = P-K 

Figure 4. Self-energy diagram involving quarks loop in HTL-approximation. 

4.2 Hadronic spectral and correlation function in HTL-approximation 

The hadronic spectral functions, <7m(w,p) in (34), of the temporal correlators are pro- 
portional to the imaginary part of the quark loop diagram. The meson spectral func- 
tion [21], constructed from two quark propagators, will have pole-pole, pole-cut and 
cut-cut contributions, <t m (w) = cr pv (uj) + a pc (uj) + cf cc (uj). Free meson spectral 
functions are also obtained using bare propagators and vertices in figure 4 as [21,22] 
a^ c (lj) = j^u 2 tanh (^) om where om = (+) — 1 for (pseudo)scalar and (— ) + 2 
for (pseudo)vector. 




012345 01234 5 

Figure 5. (a) Pseudoscalar and (b) vector meson spectral functions in HTL approxi- 
mation for m q /T = 1. The free one is represented by crosses. 

In figure 5 we show the different contributions for pseudoscalar and vector mesons [21] 
for m q /T = 1, i.e., g = \/6- Generically the pole-pole term of the meson spectral function 
in HTL-approximation describes three physical processes: (1) the annihilation of collective 
quarks, (2) the annihilation of two plasminos, and (3) the transition from upper to lower 
branch. The transition process starts at zero energy and continues until the maximum 
difference a; = 0A7m q between the two branches at k — 1.18m g . At this point a Van Hove 
singularity is encountered due to a diverging density of states for the third process as noted 
above. The plasmino annihilation starts at u = 1.86 m q with another Van Hove singularity 
corresponding to the minimum of the plasmino branch at k = 0.41 m q , where again the 
density of states diverges corresponding to second process. This contribution falls off 
rapidly due to the exponentially suppressed spectral strength of the plasmino mode for large 
energies, where only the first process, quark-antiquark annihilation starting at u = 2m q , 
contributes. For large energies this dominates and approaches the free results (crosses) for 
io » m q . The pole-cut and cut-cut contributions, which involve external gluons as can 
be seen by cutting the HTL quark self energy in figure 4, lead to a smooth contribution 
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to the spectral function. The pole-pole and pole-cut contributions in both channels are of 
similar magnitude, while the cut-cut contribution at small u vanishes in the pseudoscalar 
channel but diverges in the vector channel. This has its origin in the structure of HTL 
quark-meson vertex [19,11], which is required for the vector channel containing a collinear 
singularity, whereas bare vertices are sufficient for the (pseudo)scalar channel. In a very 
recent calculation [23] using next to leading order HTL approximation, the pseudoscalar 
spectral function is found to remain almost same as the leading order results. 



4.3 Dilepton production rate in HTL- approximation 

The static dilepton production rate (p = 0) [19] is related to the vector meson channel as 

a v (u) = -ImxvH = (e*- - l)u?-^—(p = 0). (36) 

7r oa <rxdujd p 



\ LGT T-1.5T C 


T 3 0T . - 


Born, m T C O 




,,>X:--- m/T=1.0 , 

/X^N. HTL 








co/T "' 



2 4 6 8 10 



Figure 6. The dilepton production rate in HTL-approximation for m q /T — 2 is com- 
pared with Lattice data. 



In figure 6 we show the dilepton rates calculated from (36) in HTL and also in Born 
approximations, which are compared with the first dilepton rate obtained in lattice [24]. 
The comparison of lattice rate with HTL and Born rate shows that for all energies u/T > 4 
the difference is very small. For energies u/T < 3, the lattice dilepton rate drops rapidly 
and reflects sharp cut-off found in reconstructed spectral function [24] and differs from that 
of HTL rate. 



4.4 Meson correlation functions in HTL-approximation 

The temporal correlators for pseudoscalar and vector mesons can be obtained from their 
spectral functions according to (35) in order to take into account of medium effects from 
the quark propagator. At larger temperatures it seems that the pseudoscalar correlator only 
slowly approaches the free correlation function and still differs from the lattice results [24]. 
In other words, the HTL medium effects may not be sufficient to explain the deviation from 
free one. Details of the vector meson correlation function can be found in Ref. [21]. 
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Figure 7. The pseudoscalar meson correlation functions in HTL-approximation (left) 
for m q /T — 1 and compared with lattice results. 



5. Susceptibilities 



Dynamical properties of a many particle system can be investigated by employing an exter- 
nal probe, which disturbs the system only slightly in its equilibrium state, and by measuring 
the response of the system to this external perturbation. Usually one probes the dynamical 
behavior of the spontaneous fluctuations in the equilibrium state. In general, spontaneous 
fluctuations are related to correlation functions. They reflect the symmetries of the system, 
which provide important inputs for quantitative calculations of complicated many-body 
systems. Recently, screening and fluctuations of conserved quantities have been proposed 
as a probe of the quark-gluon plasma (QGP) formation in ultrarelativistic heavy-ion colli- 
sions. 



5.1 Definition 



Let O a be a Heisenberg operator. In a static and uniform external field T a , the (induced) 
expectation value of the operator O a (0, it) is written as 

TrfO (0 it) e -/3(«+«e*)l i r 
K = (O a (0, *)> F = ^pfe^j J = \ J (O a (0, *)) , (37) 

where the translational invariance is assumed and TC ex is given by TL ex — 
— ^ Q J d 3 x O a (0, it) T a . The (static) susceptibility \afi is defined as 



= P [ d 3 x (O a (0, It) O a (0, if) ) , (38) 



assuming no broken symmetry (O a (0, it)) = (O a (0, Tj )) = 0. (O a (0, x)O a (0, 0)> is 
the two point correlation function with operators evaluated at equal times. 



5.2 Quark number susceptibility: 



The quark number susceptibility is the measure of the response of the quark number density 
with infinitesimal changes in the quark chemical potential [i q + 8[i q and can be written as 
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^ - 1: 



(3 



J d 3 x (jo(0, f)io(0, 0)> = p J d 3 x S o(0, x) , (39) 



where 5*00(0,5?) is the time-time component of vector meson correlator S^ u (t,x) = 
{jfi(t, x)j u (0, 0)). With the Fourier transform of Soo(0, x) the (39) becomes 

— 5 00 (w,0) with S 00 (u;,p) = — — ^Imn 00 (w,p) . (40) 

Using the imaginary part of the 1-loop vector meson self-energy within the HTL- 
approximation given in figure 8 the quark number susceptibility [25] is computed through 
(40). It is found to contain [25] only the quasiparticle (QP) contributions following from 
the QP dispersion relation D± = 0, given in (30). Note that the transversality condition in 
figure 8 eliminates the Landau damping (LD) contributions originating from the space like 
region. 





Q=P-K 



Figure 8. The 1-loop vector meson self energy diagram. 




°" 6 1 1.5 2 2.5 3 3.5 4 

Figure 9. The quark number susceptibility in HTL approximation scaled with free 
susceptibility is plotted as a function of T/Tc with Tc/^ms = 0-69. In each band 
the lower curve corresponds to the choice of the renormalization scale ft — 2nT and 
the upper one to p, — 2-kT. The squares represent recent lattice data. 

The quark number susceptibility is plotted in figure 9 and compared with the lattice re- 
sults. The result (band between the solid line denoted by QP) agrees well with recent lattice 
data [26]. (Here we used the 2-loop result for the running coupling constant g(T). For the 
renormalization scale p, entering the running coupling constant, we choose two different 
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values leading to the bands.) However, it overcounts the leading order perturbative result 
(dotted band denoted by PLO) given as Xq/Xf = (1 — 2a s /7r + • • •). We note that this can 
be cured by going to the calculation at 2-loop order. Recently, it has been shown that the 
2-loop order calculation [18] for the thermodynamic potential reveals the correct inclusion 
of the leading order effects but happens to be very close to the 1-loop result. Now, em- 
ploying only resummed propagators and bare vertices, the QNS contains LD contributions 
coming from the discontinuity of the 2-point HTL function in addition to the QP contri- 
bution. The result is shown in figure 9 by QP+LD. There are also calculations for quark 
number susceptibility within the approximately self-consistent HTL approximations [27]. 



5.3 Chiral susceptibility 

The chiral susceptibility measures the response of the chiral condensate to the infinitesimal 
change of the current quark mass m + dm. The static chiral susceptibility can be obtained 

as 

//■-(-OO l 
d 3 xS(0,x) = f3 ^S(w,0), (41) 

where S(w, 0) is the Fourier transformed correlator of the scalar channel. 

Note that, in contrast to the quark number susceptibility [25] where the zero temperature 
contribution vanishes due to the transversality properties of the vector channel, the static 
chiral susceptibility contains a quadratic ultraviolet divergence coming from the zero tem- 
perature contribution. Using dimensional regularization the temperature independent term 
disappears, as there is no scale associated with it, leading to xt (T) = — Nf ^ c T 2 , where 
Nf is the number of two light quark flavours and Nq is the number of colour. The next 
order contribution follows from the two-loop thermodynamic potential given in Ref. [4]. 
However, the second derivative of this expression with respect to m diverges at to = 0. 
Hence, the static chiral susceptibility cannot be calculated consistently in usual perturba- 
tion theory beyond leading order, but requires HTL resummation [28]. 



Xc(T) 



d(qq) 
dm 



Figure 10. The 1-loop scalar meson tadpole diagram. 

Now, we compute the chiral condensate (qq) from the tadpole diagram in figure 10, 
where we use the effective HTL quark propagator, yielding [28] 
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Note that the m 2 q term is the first finite a s correction to the free chiral susceptibility, since 
the two-loop correction in usual perturbative QCD diverges for m = 0. In figure 11, we 
show the chiral susceptibility in HTL approximation. In lattice simulations a peak around 
the critical temperature associated with chiral symmetry restoration has been observed. In 
contrast to lattice QCD the HTL approximation, although taking into account the medium 
effects in plasma due to the interactions, does not contain any physics related to the chiral 
restoration, which is truly a nonperturbative effect. However, we observe a strong increase 
of the magnitude of HTL chiral susceptibility towards low temperature similar as in lattice 
simulations above Tq- 



X 
I 




x 
I 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

T (GeV) T (GeV) 

Figure 11. The HTL static susceptibility (free susceptibility subtracted) as a function 
of temperature for A-^g = 300 MeV, and Nf — 2 with the choices of the renormaliza- 
tion scale fi = 2-kT and 47rT for up to order m 2 q (left) and up to order m\ (right). 



6. Dynamical Screening 

Screening of charges in a plasma is one of the most important collective effects in plasma 
physics. In the classical limit in an isotropic and homogeneous plasma the screening po- 
tential of a point-like test charge Q at rest can be derived from the linearized Poisson 
equation, resulting in Debye screening. Earlier in most calculations of the screening po- 
tential in the QGP, the test charge was assumed to be at rest. However, quarks and gluons 
coming from initial hard processes receive a transverse momentum which causes them to 
propagate through the QGP. In addition, hydrodynamical models predict a radial outward 
flow in the fireball. Hence, it is of great interest to estimate the screening potential of a 
parton moving [29] relatively to the QGP. Chu and Matsui [30] have used the Vlasov equa- 
tion to investigate dynamic Debye screening for a heavy quark-antiquark pair traversing a 
quark-gluon plasma. They found that the screening potential becomes strongly anisotropic. 

The screening potential of a moving charge Q with velocity v follows from the linearized 
Vlasov and Poisson equations as 
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exp [— ik ■ (r — wt)] 
/c 2 Rc[e;(cj = k ■ v, k)] 



(43) 



The dielectric function following from the semi-classical Vlasov equation describing a 
collisionless plasma is related to the high-temperature limit of the polarization tensor. For 
example, the longitudinal dielectric function following from the Vlasov equation is given 
by 



n 00 (w,fc) /" 

€l((J,k) = 1 p — = 1 + 



2 

_D 

/c 2 



1 



—— In — 

2fc w 



fe 



(44) 



where the only non-classical inputs are Fermi and Bose distributions instead of the Boltz- 
mann distribution. The gluon self-energies derived within the hard thermal loop approxi- 
mation have been shown to be gauge invariant and the dielectric functions obtained from 
these are therefore also gauge invariant. 

In the general case, for parton velocities v between and 1, we have to solve (43) 
together with (44) numerically. Since the potential is not isotropic anymore due to the 
velocity vector v, we will restrict ourselves only to two cases, r parallel to v and f perpen- 
dicular to v, i.e., for illustration we consider the screening potential only in the direction 
of the moving parton. The details of the perpendicular case can be found in Ref. [29]. 




(r-vt) [fm] 

Figure 12. Screening potential parallel to the velocity of the moving parton in a QGP. 



In figure 12 the screening potential <f>/Q in w-direction is shown as a function of r' = 
r — vt, where r — \r\, between and 6 fm for various velocities. For illustration we have 
chosen a strong fine structure constant a s = g 2 / (4n) = 0.3, a temperature T — 0.25 GeV, 
and the number of quark flavors Nf = 2. The shifted potentials [28] depend only on v and 
not on t as it should be the case in a homogeneous and isotropic plasma. For r' < 1 fm one 
observes that the fall-off of the potential is stronger than for a parton at rest. The reason 
for this behavior is the fact that there is a stronger screening in the direction of the moving 
parton due to an enhancement of the particle density in the rest frame of the moving parton. 

In addition, a minimum in the screening potential at r' > 1 fm shows up due to the 
loss of spherical symmetry of the Debye screening cloud. For example, for v — 0.8 this 
minimum is at about 1.5 fm with a depth of about 8 MeV. A minimum in the screening po- 
tential is also known from non-relativistic, complex plasmas, where an attractive potential 
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even between equal charges can be found if the finite extension of the charges is considered 
[31]. A similar screening potential was found for a color charge at rest in Ref. [32], where 
a polarization tensor beyond the high-temperature limit was used. However, this approach 
has its limitation as a gauge dependent and incomplete (within the order of the coupling 
constant) approximation for the polarization tensor was used. Obviously, a minimum in 
the interparticle potential in a relativistic or non-relativistic plasma is a general feature if 
one goes beyond the Debye-Hiickel approximation by either taking quantum effects, finite 
velocities, or finite sizes of the particles into account. 

The modification of the confinement potential below the critical temperature into a 
Yukawa potential above the critical temperature might have important consequences for 
the discovery of the QGP in relativistic heavy-ion collisions. Bound states of heavy quarks, 
in particular the J/ip meson, which are produced in the initial hard scattering processes of 
the collision, will be dissociated in the QGP due to screening of the quark potential and 
break-up by energetic gluons. On the other hand, the formation of colored bound states, 
e.g., qq, qq, gg, of partons at rest has also been claimed [33] above the critical temperature 
(2T C - 3T C ) by analyzing lattice data, indicating that the plasma behaves as a strongly cou- 
pled quark-gluon plasma. The dynamical screening will have important consequences on 
the J/ty dissociation and other binary states formed in the QGP just above Tc- 

7. Summary 

In this lecture the thermal field theory within the imaginary time formalism has been in- 
troduced and some of the problems in computing physical quantities in bare perturbation 
theory are pointed out. In order to overcome these problems, the prescription of HTL re- 
summation technique is discussed and based on it an improved perturbation theory (HTLpt) 
is also briefly outlined. Using this improved HTLpt various physical quantities, hadronic 
spectral and correlation functions, quark number and chiral susceptibilities, and the dy- 
namical screening potentials have been computed and discussed at length. 
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